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Abstract 

We consider the conditions under which a rotating magnetic object can produce a magnetically powered outflow in an initially 
unmagnetized medium stratified under gravity. 3D MHD simulations are presented in which the footpoints of localized, arcade- 
shaped magnetic fields are put into rotation. It is shown how the eff'ectiveness in producing a collimated magnetically powered outflow 
depends on the rotation rate, the strength and the geometry of the field. The flows produced by uniformly rotating, non-axisymmetric 
fields are found to consist mainly of buoyant plumes heated by dissipation of rotational energy. Collimated magnetically powered 
flows are formed if the field and the rotating surface are arranged such that a toroidal magnetic field is produced. This requires 
a diff'erential rotation of the arcades' footpoints. Such jets are well-collimated; we follow their propagation through the stratified 
atmosphere over 100 times the source size. The magnetic field is tightly wound and its propagation is dominated by the development of 
non-axisymmetric instabilities. We observe a Poynting flux conversion efficiency of over 75% in the longest simulations. Applications 
to the collapsar model and protostellar jets are discussed. 

Key words. Magnetohydrodynamics (MHD) - ISM: jets and outflows - ISM: Herbig-Haro objects - Galaxies: jets - Gamma rays: 
bursts 



1. Introduction 

Models and simulations of jets produced by rotating magnetic 
fields generally assume the existence of an ordered, axially sym- 
metric large-scale field of uniform polarity anchored in the cen- 
tral engine, starting with the original models of jets from ac- 
cretion disks by | Bisnovatyi-Kogan & Ruzmaikin| (|1976I) and 
[Blandford (1976 ), or the magnetic supernova model of LeBlanc] 
[^Wils on (1970). While such ordered fields are the most eflTec- 
tive in producing jets, the question whether they actually exist in 
accretion disks or the core of a star is still quite open. The largest 
scale at which magnetorotational (MRI) turbulence in accretion 
disks shapes the magnetic field structure is set by the disk thick- 
ness, which in turn is much smaller than jets. Collapsar cores 
are also small compared to the expected GRB jet, and it is not 
at all clear why the magnetic fields there should be ordered and 
axisymmetri c. Large-scale fields are not easily trapped by an ac- 
cretion disk ( jvan Ballegooijen|1989| ). Without a large-scale field 
protruding from the disk, one may still hope to launch outflows 
by twisting fields inside the disk or by magnetic loops that ex- 
tend into the di sk corona as e.g. in|Galeev et al.| ( [1979| ); |Tout & 
|Pringle| ( [T996l ); |Uzdensky & Goodman] ( |2008| ). 



the largest visible structures, but still considerably larger than 
the central engine. 

There is some numerical evidence that small-scale fields 
can also be used to generate jets. Most of these simula- 
tions are promising in terms of outflow production but lim- 
ited to the immediate surroundings of the outflow-forming disk. 
Axisymmetric simulations of outflows generated with small 
magnetic loops were done by [Romanova et alT (1998]); 



Turner 



The huge range of involved length scales is a major issue in 
jet modeling, as numerical simulations that cover all scales are 
not feasible to date. Protostellar jets may be several parsecs long, 
launched by disks with sizes of -100 AU, i.e. about a factor 10"^ 
smaller (e.g. Shepherd et al. 2001 ). Assuming that the "launch- 
ing scale" is much smaller than the disk, perhaps on the order 
of ~1 AU, one obtains an even bigger contrast in length scales. 
Supposing that jets in AGN are launched at a few Schwarzschild 
radii from the central black hole and taking Cygnus A as an ex- 
ample, one obtains a ratio of 10^ between the jet length an d the 
size of the engine ( Krichbaum et al . 1998 ; Tadhunter et al. 2003 ). 
The final jet properties are determined before it becomes ballis- 
tic, probably at scales which are somewhat smaller than that of 



[eTaLl ( p^99l ); |Kudoh et al.1 d2602|). In their 3D simulations of 
accretion flows, |Kato et al. ( |2004| ) demonstrated that an ini- 
tially poloidal magnetic field confined within a rotating torus 
surrounding the accreting black hole can give rise to a transient 
outflow driven by accumulated toroidal fields in the form of a 
"magnetic tower" ( |Lynden-Bell|2003| ). pe Villiers et al.|p005| ) 
showed in simulations that loops of poloidal field in an accret- 
ing torus may give rise to a large-scale poloidal field as the field 
lines are stretched out in an axial outflow. The inflation and dis- 
ruption of a magnetic loop outside a disk, caused by the gen- 
eration of a toroidal field through diflferential rotation, was ob- 
served by Fendt ( 2009 ) in simulations of outflows from star-disk 
magnetospheres. The generation of magnetic flows within un- 
magnetize d surroundings has also been studied in laboratory ex- 
periments ( Hsu & Bellan|2005l Ciardi et al.|2009| ); such jets are 
strongly afl'ected by current-driven instabilities. 

The aim of the calculations presented here is to see whether 
small-scale fields in the form of loops anchored in a rotating disk 
can be used to produce jets of significant length (compared to 
the size of the source). In the cases studied the flows propagate 
and are confined in an external unmagnetized atmosphere (as op- 
posed to MoU et al.,200 8 and Moll^2009 , hereafter Papers I&II, 
where we studied jets embedded in a large scale magnetically 
dominated environment). While still idealized, this addresses en- 
vironments like protostellar jets launched into a dense cloud, 
or GRB jets launched by a collapsar core. To be investigated 
here are the circumstances under which jetlike flows are formed 
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Figure 1. Considered magnetic field geometries. In case (a), the 
field lines are one-sidedly anchored in a rigidly rotating disk. 
In case (b), the field lines' footpoints are either both anchored 
in the disk, or they are both anchored outside of it. In case (c), 
the field emerges from a limited area inside the disk which is 
diflTerentially rotating. 



that penetrate through the atmosphere instead of dissipating in it. 
As this is also a question of (non-axisymmetric) stability, three- 
dimensional simulations are necessary. One of the questions to 
be answered is whether models of jets from small-scale mag- 
netic fields are a viable alternative to those based on the twist of 
large-scale fields. 



2. Models 

The primary ingredients in all our models are a rotating disk, 
implemented as a boundary condition, and magnetic field loops 
anchored in and sticking out of the disk. The magnetic field ge- 
ometries considered are sketched in Fig.[T] In case (a), some of 
the loops have one footpoint inside the disk while the other is 
anchored outside. The resulting shear motion of the footpoints 
generates a toroidal magnetic field. In case (b), the loops are ar- 
ranged such that they are not sheared. Here, a toroidal field can 
only be produced by the inertia of the material above the disk, 
such as in conventional models of jets from large-scale magnetic 
fields. This case would be a model for a (solidly) rotating stellar 
core inside a nonmagnetic envelope. Finally, we consider a case 
in which all loops are anchored in the disk and shear is created 
through differential rotation (c). This case would be more repre- 
sentative of small-scale fields generated in an accretion disk. 

The magnetic field loops are established either by a suitable 
potential field, imposed as initial condition (setup D), or by con- 
tinuous "injection" of the field from below the disk (setup E), 
see Fig. |3] and Sect. |3.1.2| for details. While the magnetic field 
geometry is similar in both cases, setup E (for "Emerging field") 
is meant as an idealization of the field loops emerging from mag- 
netic turbulence in an accretion disk. 



The model chosen as initial condition for the atmosphere is 
a hydrostatic equilibrium stratification in the gravitational po- 
tential of a point mass (the origin). Temperature T and density 
p vary with distance r from the point mass as r"^ and r"^, re- 
spectively. The temperature is thus a constant fraction of the 
virial temperature, as in advection-dominated accretion flows. 
The density profile has the property that the amount of mass in a 
cone centered on the origin is constant per decade in distance r. 
The stratification is thus scale-free, and a prospective jet encoun- 
ters the same amount of atmosphere mass per decade traveled. 
With the equation of state used in the calculations, an ideal gas 
with ratio of specific heats y = 5/3, the stratification is convec- 
tively stable. 




Figure 2. Sketch of the coordinate system and the computational 
grid. The jets propagate in equatorial direction of the grid coor- 
dinates (r, 6, (p). To describe the results, the alternate spherical 
coordinate system (r, (f) is used, which takes the jet's central 
axis as the polar axis. 



3. Methods 

3.1. Numerical MHD solver, grid and coordinates 

We numerically solved the ideal adiabatic MHD equations with 
a static external gravitational potential O oc r"^ on a spherical 
grid (r, 6, cp). The jets propagate in equatorial direction along the 
y-Sixis (6 = (p = 7r/2), about which the computational volume 
covers a range = Acp in the angular directions, see Fig. [2] The 
spacing is uniform in 6 and and logarithmic in r. Such a grid 
is more economical than a Cartesian one since it expands in jet 
direction. Unlike a spherical grid in polar direction z, it is free of 
singularities which may cause numerical problems or artifacts. 

As jet physics is better described in a coordinate system 
which takes the jet's axis y as the polar axis, we introduce an- 
other spherical coordinate system (r, i^, (f) for that purpose, see 
also Fig. |2] /? := r sin ?^ denotes the orthogonal distance to that 
axis (cylindrical radius). 

The simulations were performed with a newly developed 
Eulerian MHD code (Obergaulinger'2008 ). It is based on a flux- 
conservative finite-volume formulation of the MHD equations 
and the constraint tran sport scheme to maintain a divergence- 

I [1988]) 



free magnetic field (Evans & Hawley [1988 ). Using high- 
resolution shock capturing methods (e.g., LeVeque|1992 ), it al- 
lows a choice of various optional high-order reconstruction al- 
gorithms and approximate Riemann solvers based on the multi- 
stage method (Tor o & Titarev|2006| ). The simulations presented 
here were performed with a fifth order monotonicity -preserving 
reconstruction scheme ( [Suresh & Huynh|[l997| ), together with 
the HLL Riemann solver ( Harten|1983| ) and third order Runge- 
Kutta time stepping. 

As is done in other astrophysical MHD simulations, the im- 
plicit assumption has been made that the magnetic reconnection 
in the settings studied is of the "fast" type, for which the recon- 
nection implicit in the numerical diffusion of the code may be a 
fair representation of the subgrid dynamics. 



3.1 .1 . Setup with an arcade-shaped, exponentially 
decreasing initial field (D) 

The initial condition for the magnetic field in this case is a 
two-dimensional potential field which is independent of the x- 
coordinate (see Figs. |2|3| for the coordinate system) and falls oflT 
exponentially with height It is generated by the vector poten- 
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Table 1. Normalization units 



Figure 3. Schematic of simulation setup E described in 
Sect. 3.1.2[ in the variant sketched in Fig.[T^. A vertical velocity 
field (green arrows) pushes the magnetic field (blue lines) from 
the lower boundary (r = rb) into the computational volume. An 
azimuthal velocity field within R < R\y twists the created mag- 
netic arcades (dotted blue lines) which are one-sidedly anchored 
in the rotating region. 



tial 



A = ^b^sinQjexp^^^^-^j^;,, 



(1) 



where An is the width of the arcades and A is the scale height of 
B = \B\. is the field strength at the intersection of the lower 
boundary r = rb with the central axis (y). To avoid numerical 
problems at the lateral boundaries, we limit the number of arcs 
to 4, setting A = B = Ofor\z\> 2A7[ and for \x\ > 2An. 

The magnetic field is embedded in a spherically stratified 
atmosphere with p oc and p oc r~^, held in hydrostatic equi- 
librium by the static gravitational field O oc r"^ of a point mass 
at the coordinate origin. Temperature and sound speed vary as 
T oc and Cs oc r"^/^, respectively. In the outer, unmagnetized 
region (\x\ , |z| > 2 An), we compensated for the absence of mag- 
netic pressure by increasing the gas pressure. To maintain hydro- 
static equilibrium, the density is also increased correspondingly. 

At the lower (r = rb) boundary we maintain, through "ghost 
cells" outside of the computational domain, an azimuthal veloc- 
ity field V = v^e^p corresponding to rigid rotation: v<^ = v^^R/R\) 
for R < Rh and elsewhere. All quantities except for B are fixed 
at their initial values in the ghost cells; B is extrapolated from 
the interior of the domain. At the sides (6 and 0) and top (upper 
r) of the domain, we use open boundary conditions which allow 
for an almost force-free outflow of material and cause no evident 
artifacts in the form of reflections. 

A is chosen such that in the x = plane By is positive for 
\z\ < Rh and negative for 7?b < \z\ < 27?b, i.e. arcades which 
start at \z\ < Rh from the r = rb surface have their second foot- 
point outside the rotating region. This is achieved by choosing 
the length scale of the arcade as /I = 27?b /n. 

3.1 .2. Setup with magnetic field arcades emerging from the 
bottom boundary (E) 

In this case, the initial condition is again an equilibrium stratifi- 
cation with p oc , p oc and O oc r"^ However, unlike in 
setup D, the atmosphere is completely unmagnetized. 

The magnetic field enters the domain through the lower 
boundary r = r^, the conditions of which are determined through 
ghost cells at r < rb. There, we impose a transverse magnetic 
field Bq Bz) with constant amplitude ^b and a polarity that 
alternates with z = rb cos 6 in step sizes of 27?b, the diameter of 
the rotating surface (described below). Where the polarity of Bq 
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changes, the solenoidality of the magnetic field is maintained by 
an appropriate Br. The field lines thus have the shape sketched in 
Fig. [3] An approximate equilibrium is maintained through low- 
ering the gas pressure by the value of the magnetic pressure, as 
far as this is possible (p must not be negative). For this to work, 
y^b •= ^^Pb/^l must be greater than one, which limits the pos- 
sible field strengths to ^b < 1 in our system of units (described 
below). 

To model the emergence of magnetic fields into the atmo- 
sphere we impose the radial velocity field 
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in the ghost cells of the lower boundary. The injection velocity is 
too small to form a jet by itself: the maximum amounts to 41% 
of the escape velocity at the boundary, gas with this speed gets 
theoretically as far as r 1.2rb without acceleration other than 
gravity. In addition to the radial velocity field, we maintain an 
azimuthal velocity field v^p oc R (rigid rotation) within R < R\y. 
The shape of the emerging magnetic field is sketched in Fig.[T^. 
For the other boundaries, we use the same outflow conditions as 
in the model described in the preceding section. 

Two variations of the above-described setup have also been 
studied, corresponding to the cases (b) and (c) described in 
Sect. [2] In the first, the positions where Bq changes its polar- 
ity, and with it the radial velocity field that injects the magnetic 
field, are shifted by 7?b in -z direction. Thus, both footpoints of 
an individual field loop rotate with the same angular velocity. 
In the second case (c), the magnetic field and vertical veloc- 
ity in the lower boundary are changed such that only a small, 
confined arcade in 0.2 < R/Rh < 0.8 with a width of 0.27?b 
emerges, and the azimuthal velocity field has a Keplerian profile 
with = v^l^ iOAR^/R for 0.1 < R/R^ < 1. The inner edge 
(R = 0.27?b) of the arcade rotates twice as fast as the outer edge 
(R = 0.87?b), the diflTerence in rotation velocity is v^^/ VS. 



3.2. Parameters and units 

The models described above contain the following 6 parameters, 
not all of which are independent: rb, 7?b, Pb, /^b, and v^^^. 
We eliminate the dependences by taking /q = 27?b, Po = Pb and 
Pq = p^ as units of length, density and pressure, and expressing 
all physical quantities in terms of these, see Table[T] The remain- 
ing 3 parameters can then be expressed as dimensionless num- 
bers. The first of these is /o/^b, which is a measure for the cur- 
vature of the rotating surface. In all simulations presented here, 
this curvature is small and probably does not influence the results 
significantly. The parameter has technical significance, however, 
because it also controls the opening of the lateral boundaries. 
The remaining two numbers are chosen to be y^b •= ^^Po/B^, 
which controls the strength of the magnetic field, and the Mach 
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Figure 4. Selected magnetic field lines in a simulation with an 
arcade-shaped initial field (setup D). The lines are closing with 
the lower boundary, ascending inside the jet and descending at 
the outside. 




Figure 5. Left: Magnetic field and temperature in a horizontal 
slice next to the bottom boundary in a simulation with an arcade- 
shaped initial field (setup D) with shear (Fig. [T^) after 1 (bot- 
tom left panel) and after 32 revolutions {top left panel). The tex- 
ture shows the geometry of the parallel field component and the 
green/yellow coloring gives the magnitude of the normal com- 
ponent. Regions of increased temperature are tinted in red (low 
temperatures are transparent). Right: Evolution of the radial ve- 
locity in a meridional slice. The jet is eventually "choked oflT" at 
the bottom. 



number Mb := v^^^/cso, which determines the speed of rotation. 
Note that in the emerging field model (E), is a measure of the 
field strength but not a local plasma-beta value, since the pres- 
sure entering its definition is not measured at the same location 
as the field strength. For the sake of clarity, we usually omit the 
units in the presentation of the results. The concerned quantity 
is then measured in terms of the associated normalization unit 
listed in Table [T] 



4. Results 

As explained in Sect. [2] above, the calculations were done with 
diflTerent models for the rotating magnetic field configuration. It 
turns out that some of these cases produce long-lived jets, oth- 
ers only transient flows or flows that dissipate in the atmosphere 
close to the source. The results presented i n the f ollowing sub- 
section were obtained with setup D (Sect. |3.1.1| ), those in the 
subsequent subsections were obtained with setup E (Sect.|3.1.2|). 



4.1. Transient jets witli arcade-sliaped initial field 

Simulation setup D with the field geometry in Fig.[T^ was found 
to produce extended, collimated outflows. However, the mag- 
netic field at the base of the jet decays in this case (hence the 
"D") and is not replenished, for which reason these jets are not 
permanent. 

The forced rotation at the lower boundary stretches the mag- 
netic field lines in azimuthal ((p) direction around the axis of ro- 



tation (y). The gas is accelerated upwards and the magnetic field 
assumes a tangled helical structure. A good momentary accel- 
eration was obtained with the parameters (J3\y = 1/4, Mb = 18). 
The magnetic field in this case is depicted in Fig. |4] The jet di- 
ameter is resolved with about 35 grid cells in this simulation. 
The jet attains a height of about 30 times its initial diameter, the 
factor of expansion being about 7. It reaches velocities that are 
close to Mb and about a factor of 10 above the escape speed. The 
field lines are closing with the lower boundary: inside the jet, 
the radial magnetic field has the same polarity as on the disk, 
whereas the net radial magnetic flux J BrdA, integrated over the 
r = const surface, is virtually zero at all times and all radii. 

The "absolute flux" J \Br\ dA decreases linearly with time 
near the lower boundary. Increased temperatures near the out- 
side of the rotating surface indicate that a substantial amount of 
magnetic field is being dissipated there, see Fig. [5] Due to the 
rotation, the magnetic field assumes a vortex-like structure, at 
the border of which magnetic field lines of opposite polarities 
become entangled. This leads to a continuous decay of the mag- 
netic field. Such a decay of the magnetic field, through wrap- 
ping up of field lines followed by cancellation through diflfu- 
sion, is called "convective expulsion" in overturning flows (e.g. 
Zerdovich|1956[ [Parker|1963] ). The magnetic field in the phys- 



ical source of a successful jet evidently must be of a diflTerent 
nature. This calls for a modification of the boundary conditions 
at the base of our simulations. Such cases are discussed in the 
following sections. 
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Figure 6. Top panel: Maximum velocity over escape speed for 
different parameters (y^b, Mb) in otherwise similar simulations of 
the emerging field model (E), variant (a). The strength of the jet 
decreases with a decrease of the field strength (increasing y^b, 
represented by different line styles) and with a decrease of the 
rotation velocity (Mb, represented by different colors). Bottom 
panel: Strength of the jets as a function of the rotation rate and 
field strength parameters. 



4.2. Jets from emerging fields: parameter study 

Simulation setup E, in which the magnetic field loops above the 
rotating surface are continuously replenished, produced long- 
lasting jets that propagate considerable distances. Before pre- 
senting big simulations with long jets in the next two sections, 
we discuss the influence of the parameters (J3\y, Mb) on the re- 
sults by comparing a series of smaller, computationally cheaper 
simulations in which these parameters are varied. Although the 
proximity of the upper boundary possibly influences the results 
(see discussion of instabilities below), the small simulations give 
clear indications about which parameters yield efliicient jets. 

The simulations cover the relatively short distance of 20 
length units (40 disk radii), with 10 < r < 30 and 6,(p = 7i/2 ± 
71 19, the resolution being 256 x 96 x 96. The atmospheric den- 
sity differs by a factor 27 between the lower and upper bound- 
aries. The simulations were pursued until t ^ 100. In the cases 
where a jet was successfully launched, this sufl&ces for about 2-3 
passages through the computational volume. The computational 
cost for this was 1-3 wallclock days on 36 processors with MPI 
parallelization. 

To test the flow's ability to penetrate the atmosphere in the 
setup with shear, Fig.[T^, we compare the jet velocities halfway 
through the simulated distance against the escape velocity in 




Figure?. Parallel magnetic field (wedges) and temperature 
(color) in a meridional slice through the jet in a simulation with 
shear (top panel) and in two simulations without shear {mid- 
dle/bottom panel). Collimated, magnetically driven jets form 
only in cases where a toroidal field develops through a shear of 
the magnetic field. 



Fig. [6] The jet velocity decreases with both decreasing rotation 
velocity Mb and decreasing magnetic field strength ^b (or in- 
creasing y^b)- It is, however, much more sensitive to the former 
parameter: y^b must be changed by orders of magnitude to get a 
significant impact on jet velocity whereas with Mb, a factor of or- 
der unity sufl&ces. The bottom panel in Fig.[6]shows which flows 
pass this test, with those exceeding the escape speed marked 
in green, the others in red. Failed jets do not reach the upper 
boundary. The background atmosphere, whose equilibrium is 
perturbed, tends to fall down on them. The fixed conditions at 
the lower boundary avert a pile-up of thermal energy: downward 
flowing hot gas vanishes across the boundary and the injected 
gas has a constant temperature. A delayed onset of a flow due to 
accumulated heat is unlikely for this reason. 

There are instabilities in all cases where a jet develops. 
The instabilities develop mainly after the first passage through 
the computational domain, in the form of helical displacements 
and/or a change of direction of the whole jet by up to several de- 
grees. The latter are likely modes with wavelengths longer than 
the computational domain, i.e. one sees only the lower part of 
what would be a kink if the radial extent of the domain was 
larger. Jets created by stronger magnetic fields (J3\y < 4.9) tend 
to show this kind of incipient instability. Jets from weaker fields 
move slower and develop pronounced helical deformations al- 
ready within the computational domain. The proximity of the 
upper boundary does not allow for more conclusive statements 
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Figure 9. Meridional slices through the jet presented in Sect. 4.3 
The magnetic jet is accelerated into a non-magnetic atmosphere 
penetrating through denser material. The hot, magnetic material 
in the environment of the upper part of the jet are remnants of 
previous instabilities. 



about differences in instability behavior within the limits of this 
parameter study. A series of large and expensive simulations of 
the kind presented in the next section would be needed for that. 

If the configuration is shifted such that no shear in the mag- 
netic field occurs (Fig. [T])), a jet does not develop. The middle 
panel in Fig. |7] shows such a case. The magnetic field lines are 
unconnected to the surrounding atmosphere and the field is not 
amplified by shear. This appears to be suflftcient to allow them to 
rotate without producing a magnetically powered flow. However, 
for very large values of Mb and y^b, an outflow forms at the edge 
of the rotating disk, see bottom panel in Fig. [7] This outflow 
appears to be driven by thermal buoyancy associated with the 
dissipation of rotational energy at the disk's edge. Having the 
shape of a hollow cone with a large opening angle, it is clearly 
distinctive from the jets discussed above and the rest of this pa- 
per. 

4.3. Jets from emerging fields: a large simulation 

The long-range behavior of jets from emerging fields (setup E) 
in the case with shear (Fig. [T^) was investigated in a large sim- 
ulation that covers the radial range 10 < r < 110. Along this 
range, the atmospheric density decreases by a factor 1300. The 
parameters used were (J3\y = 4, Mb = 3). The jet crosses the up- 
per boundary at the physical time t ^ 110, which corresponds to 
105 complete rotations of the disk; the simulation was stopped at 
t = 539. For this, 19 wallclock days on 64 processors with MPI 
parallelization were needed, the resolution being 400 x 160x 160. 
The lateral boundaries are at ^, = 7r/2±7r/6; the jet stays within 
these boundaries at all times. 

The wound-up magnetic field rises upwards as shown in 
Fig. [8] with a collimated outflow forming in vertical direction. 
Not counting the changes in direction caused by instabilities, the 
jet's half opening angle is about 5°. Inside the jet, the field lines 
spiral upwards, with the radial magnetic flux having the same 




Figure 10. Spatial and time-averaged val ues o f various quanti- 

The dotted lines 



ties in the simulation presented in Sect. [43 
represent the initial atmospheric value, the 
"parallel to v". 



symbol stands for 



sign as the radial field on the rotating surface. The total flux is, in 
comparison, close to zero at most radii and times. Accordingly, 
visualizations of the field lines give the appearance that most 
lines remain connected to the lower boundary with both ends 
until the jet crosses the upper boundary. 

Fig. [To] shows the average value of various quantities across 
the flow, i.e. the average of a quantity X is calculated as 
over the r = const surface. The average ve- 
locity increases mainly below r ^ 20, at which point it exceeds 
the average fast magnetosonic speed v^^ = -h v^. The peak ve- 
locity (maximum of on r = const) is constant beyond r 30 
and is twice as large as the average velocity (4 instead of 2). The 
average density is at all radii significantly lower than the atmo- 
spheric value, the temperature (c^) is increased with respect to 
the environment. The transversal (toroidal) magnetic field dom- 
inates over the radial component at medium and large distances. 
The magnetic pressure dominates over the gas pressure by a fac- 
tor of about 2-4 at all distances in the jet. 

The dependence of the energy flow rates with distance gives 
information about energy transformations taking place in the jet. 
The total energy flow rate 6tot(^, r) is obtained from a surface 
integral over the components of the radial energy flux 
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(3) 
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Figure 8. Selected magnetic field lines in the simulation presented in Sect. |4.3[ The color gives the magnetic field energy. The 
boundary conditions create magnetic field arcades above the boundary, the arcades are twisted in the center by the imposed rotation 
and a jet with helically shaped field lines forms along the axis of rotation. The jet is subject to occasional helical deformations. 
Movies of this and other simulations can be found at http : //www . mpa- garching . mpg . de/ ~rmo/pap3/index . html 
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Figure 11. Energy flow rates in the jet presented in Sect. |4.3 
Magnetic enthalpy (integral over Poynting flux, blue line) is 
converted into kinetic energy (green line) and potential energy 
(brown line). 



(4) 



'grav 



with S r being the radial component of the Poynting vector: 

We denote the individual components of 6tot by 5kin, Sthrm, 
and 6niag in order of their appearance in Eq. ([3]). The compo- 
nents of the magnetic enthalpy flow rate are denoted by 6^^g 
in order of their appearance in Eq. (|4]). The energy flow rates 
in the present simulation are plotted in Fig. [TT] The plot shows 
how magnetic enthalpy is converted into kinetic and potential 
energy. The conversion is efl&cient in that less than 25% of the 
initial 5mag remains in the jet when it reaches the upper bound- 
ary. Thermal energy is less important, 6thrm being reduced only 
by half as many units as 6mag- The most important components 
of 6mag are 6^ag' which is associated with the advection of the 
azimuthal field and, at low radii, fi^ag' the work done by the 



azimuthal flow against the azimuthal component of magnetic 
stress. 6^ag is virtually zero. In summary, the qualitative behav- 
ior of the energy flow rates is similar to those in simulations 
of jets generated by twisting a uniformly polarized large-scale 
magnetic field (Papers I&II). 

The jet has a lot of substructure, see Fig. |9] Above r ^ 20, 
it is affected by recurrent instabilities that divert its course away 
from the j-axis by several degrees. Stirred up ambient material 
contributes to the unstable behavior. Occasionally, the jet de- 
velops a pronounced helical shape as is characteristic for kink 
instabilities. The helix makes approximately one complete turn 
within the computational volume, out of which it is rapidly ad- 
vected. For comparison, the magnetic pitch is <10 inside the jet 
(cf. Fig. [8]), which is at least a factor 10 below the wavelength of 
the prominent instability. The largest deflections from the central 
axis amount to approximately 10°. Being continually recreated 
at its base, the jet survives in a time-averaged sense. 



4.4. Jets from a magnetic arcade on a differentially rotating 
surface 

We also produced jets with emerging field loops (setup E) on a 
diflferentially rotating surface, see sketch in Fig. [T]:. These jets 
are mostly the result of thermal buoyancy, driven by dissipative 
heating of the near-disk atmosphere. Nonetheless, the presence 
of a magnetic field is required to launch a directed flow. 

The simulations cover the same physical range as the one 
presented in the preceding section. The resolution is 400 x 256 x 
256, so that the diameter of the diflferentially rotating surface is 
resolved with about 25 cells. The longest simulation, which also 
produces the most eflftcient jet, has the parameter values (J3\y = 
1.2, Mb = 15). It ran for 14 wallclock days on 96 processors 
and covers almost two jet crossings, the upper boundary being 
reached ait ^ 62. 

In the jet of the preceding section, Poynting flux was the 
main source of energy to accelerate the jet. This is clearly dif- 
ferent here. The simulation with (J3\y = 1.2, Mb = 15) attains a 
momentary quasi- stationary state, with 6tot ~ const in the en- 
ergy flow, at ^ ^ 80 ... 90. In sharp contrast to Fig. [TT] 6thrm is 
the dominating component at the beginning, being a factor of 
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Figure 12. Meridional slices through the jets presented in 
Sect. |4.4| A jet forms only if a magnetic field is present and the 
rotation velocity is sufi&ciently large. Movies of such jets can be 
found at the web address stated in Fig. [8] 



~3 larger than 6mag- The equipartition point S^^ = 6mag is very 
close to the lower boundary, the radius where 5kin = ^thrm is 
larger, at r 15. 5mag also diminishes (by about 80%), but it is 
only a small fraction of the final 6kin. 

The excess temperature responsible for the high 6thrm turns 
up near the central axis next to the lower boundary and increases 
with the speed of rotation. The heat is apparently produced by 
numerical viscosity, i.e. dissipation of kinetic energy due to the 
coarse resolution of the rapidly rotating inner part of the disk. 
Unfortunately, a much higher resolution is unfeasible and a large 
Mb is needed for an efl&cient production of B^^. The simulations 
are therefore relevant for cases in which dissipative heating is of 
importance. 

The appearance of the flow becomes more jetlike with in- 
creasing rotation velocity. In a simulation with (J3\y = 1.2, Mb = 
3), corresponding to a diflTerential rotation velocity of about 1.1 
of the outer loops' footpoints, the jet terminates at a few disk 
radii from the source in a tightly wound helix, similar to the free 
end of a garden hose. With Mb = 9, the flow has a turbulent, 
"smoke stack"-like appearance. The jet produced with Mb = 15 
is subject to non-axisymmetric perturbations, but much more co- 
herent. 

Despite the main energy source being thermal, the results 
depend greatly on whether a magnetic field is present or not. 
In a simulation with (J3\y = 10^^, Mb = 15), i.e. with virtually no 
magnetic field, there is no unidirectional flow. Rather, material is 
ejected sideways as well as in the upward direction, see Fig. 12 
The simulation soon crashes, presumably due to too much heat 
accumulating. The same simulation with /3\^ = 1.2 is well colli- 
mated and plows eflTortlessly through the dense parts of the at- 
mosphere. 



5. Summary and discussion 

We produced collimated outflows in numerical simulations by 
rotating the footpoints of arcade- shaped magnetic loops. The 



longest of these jets cross a computational domain which is two 
orders of magnitude larger than the size of the source. 

Flow acceleration via dissipative as well as non-dissipative 
processes was observed. The latter relies on the presence of shear 
in the form of diflferential rotation of the footpoints of individ- 
ual field loops, which results into the development of a toroidal 
magnetic field component. The resulting jet has a helical mag- 
netic field structure with field lines running back to the source 
outside the jet. 

The case of uniformly rotating magnetic loops has also been 
studied. It turns out to be markedly less eflfective at produc- 
ing outflows; in fact, no jetlike outflows were observed for uni- 
formly rotating arcades within the numerically accessible pa- 
rameter range. The shear between the rotating loop and the sta- 
tionary atmosphere surrounding it also produces a toroidal field 
component, but this appears to be much less eflTective than direct 
shearing of field lines inside the rotating source. 

We find that the continuation of the outflow beyond the ini- 
tial transient depends on the way the magnetic field is main- 
tained at the lower boundary. DiflTerential rotation acting on a 
non-axisymmetric field quickly dissipates the field through con- 
vective expulsion, with the result being that the jet is "choked 
oflT" at the base. We have compensated for this by adding an 
inflow of magnetic field at the base. This would represent, for 
example, the emergence of loops of magnetic field into the at- 
mosphere of an accretion disk. 

Apart from geometric requirements, the most important pa- 
rameter for efl&cient jets turned out to be the (diflTerential) speed 
of rotation of the footpoints. In the models calculated, the nec- 
essary speed is of the order of the escape velocity from the 
disk. This is more than what is needed in simulations with a 
long-scale poloidal initial field (e.g. those presented in Paper I). 
Consequently, the critical points (sonic, Alfven, fast magne- 
tosonic) and the equipartition point between the magnetic and 
kinetic energy flow rates are reached at smaller distances. The 
conversion of magnetic enthalpy (and the components thereof) 
to kinetic energy is similar apart from that, being fairly eflftcient. 

In the adiabatic calculations presented here, thermal energy 
from hydrodynamic and magnetic dissipation stays in the flow. 
This is in contrast with our earlier calculations (Papers I&II), 
where we assumed optically thin environments in which much 
of the dissipated energy is lost by radiation. The present models 
assume no energy loss and are most relevant for optically thick 
conditions. In the stratified atmosphere of these models, heat- 
ing contributes to the flow by thermal buoyancy. In some of the 
cases presented (Sect. |4.4|), this is the dominant driving mecha- 



nism. Though a wound-up magnetic field is also present in these 
buoyant plume flows, they probably do not qualify as magnet- 
ically driven. The non-axisymmetric "stirring" by the rotating 
magnetic field at the base that drives them may, however, well 
be relevant for the case of a (rapidly) rotating core inside a dense 
stellar envelope. Even though its field configuration is not of the 
right kind to produce a magnetically powered outflow, dissipa- 
tion of rotational energy by "stirring" may still have powerful 
effects on the envelope. 

The magnetic field in the jets presented here is strongly 
wound, with the field lines making many turns along the length 
of the jet, see Fig.[8]for an example. Current-driven instabilities, 
in particular kink modes, are therefore to be expected and can 
indeed be found in the simulations. The turbulent wiggling of 
the jet increases the interaction with the ambient medium, lead- 
ing to an increased entrainment of material. At large distances 
(larger than those covered by the simulations), it is possible that 
the field configuration is aflfected by instability-induced dissi- 
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pation of magnetic energy, presumably with dynamical conse- 
quences for the jet. The effect would likely be the same as in jets 
launched with large-scale magnetic fields (see Paper II for that). 
Kelvin-Helmholtz instabilities at the interface between the jet 
and the ambient medium are suppressed by the relatively coarse 
resolution used in the simulations. 

The final state of the jet, at large distances from the source, 
depends on the impact of instabilities in the region beyond the 
computational volume of the simulations presented here. It is 
reasonable to believe that the jet continues as a ballistic flow 
when the instabilities cease to be eflTective (due to the expan- 
sion of the jet and a possible decay of the toroidal field). A "dis- 
ruption" of the jet would require a strong interaction with the 
medium into which it propagates; the atmosphere at large dis- 
tances is probably too thin for that. 

Although the simulations presented here demonstrate pos- 
sible forms of rotation-induced acceleration, they probably do 
not represent very realistic models of actual accretion disks. 
One may wonder what happens in a more realistic scenario 
with more complicated or chaotically arranged magnetic field 
loops emerging from the disk. Obviously, magnetic reconnection 
would likely play a major role in such a scenario. Nonetheless, 
the surviving magnetic field is stretched out in azimuthal direc- 
tion by the rotation, and the free energy in the toroidal field is 
transformed into an outflow. Temporary fluctuations are likely 
to occur if the supply of magnetic field loops is not continuous. 
The coherent length of the jet will then be determined by the 
speed with which reconnection processes destroy the initial field 
by which the flow is produced. The jet may flare up anew when 
new loops emerge. This effect could be responsible for the tem- 
porary fluctuations in gamma ray bursts. 

Loops of magnetic field can also be found in the solar corona 
and their evolution there is likely connected with explosive ejec- 
tions of material (CMEs, see e.g. [Forbes et a l. 2006 for a review). 
One general idea is that the diflTerential rotation as well as con- 
vective motions at the solar photosphere produces shearing, and 
the energy in the sheared/twisted field is set free in an explosive 
event (e.g. Van B allegooijen & Martens||1989[ Zuccarello et al. 
2009]). In contrast to the jet models discussed here, the footpoints 
of a solar magnetic arcade is unlikely to be twisted by a (pure) 
azimuthal motion. Non-azimuthal shearing is not likely to pro- 
duce an extended helical field (as in Fig. [4]) and the plasma accel- 
eration is hence diflTerent from that in a magnetically driven jet. 
Also, it seems that the energy in CMEs is first accumulated and 
then released in a sudden event, much unlike the rather smoothly 
rising jets presented here. 
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